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Final  Performance  Report  for  Grant  FA9550- 13- 1-0033: 
“Energy  Relations  in  Natural  and  Artificial  Diamagnetic  Materials” 

Program  Manager:  Dr.  Arje  Nachman,  Principal  Investigator:  Dr.  Arthur  D.  Yaghjian 


1  Introduction 

As  a  consequence  of  the  insight  gained  from  the  rigorous  homogenization  theory  developed  by  the  principal 
investigator  in  collaboration  with  Andrea  Alu  and  Mario  Silveirinha  [1]  for  the  magnetic  permeability  as  well 
as  electric  permittivity  of  3D  spatially  dispersive  periodic  arrays,  the  AFOSR  three-year  Grant  FA9550-13- 
1-0033  for  $1 10,889  and  its  Additional  Work  Effort  for  $72,600  was  awarded  to  S4-Inc  to  have  Dr.  Yaghjian 
apply  this  theory  to  determine  consistent  expressions  for  the  power  and  energy  densities  in  a  diamagnetic 
continuum.  The  tasks  proposed  and  successfully  completed  for  the  three  years  were  in  brief: 

1 .  To  evaluate  the  external  and  internal  power  and  energy  relations  for  a  realistic  model  of  a  diamagnetic 
molecule  or  inclusion,  such  as  a  perfectly  electrically  conducting  (PEC)  wire  loop; 

2.  To  determine  a  reliable  method  for  obtaining  macroscopic  continuum  relations  of  power  and  en¬ 
ergy  from  the  microscopic  Maxwellian  equations  describing  the  power  and  energy  of  diamagnetic 
molecules  or  inclusions. 

3.  To  incorporate  the  results  for  the  macroscopic  power  and  energy  relations  into  the  continuum  formu¬ 
lation  for  diamagnetic  materials  or  metamaterials; 

4.  To  combine  the  diamagnetic  relations  for  power  and  energy  with  those  of  electric  polarization  to 
obtain  positive  semi-definite  (that  is,  nonnegative)  energy  expressions  for  diamagnetic  continua,  then 
document  the  results  of  the  research  in  archival  publications. 

These  accomplishments,  which  were  recently  submitted  to  the  journal  Physical  Review  B ,  will  be  summa¬ 
rized  in  the  remainder  of  this  final  report.  Some  of  the  results  can  also  be  found  in  the  references  [2]— [4] . 

There  are  several  reasons  for  determining  nonnegative  macroscopic  energies  from  the  microscopic 
Maxwellian  equations  for  classical  diamagnetic  dipolar  continua.  First,  it  is  one  of  the  remaining  unre¬ 
solved  problems  in  classical  electromagnetic  theory  even  though  the  problem  is  relatively  easy  to  state: 
given  a  volume  of  a  macroscopic  continuum  satisfying  Maxwell’s  dipolar  equations  and  illuminated  by  ex¬ 
ternal  fields,  are  there  macroscopic  energy  densities  that  never  become  negative  for  all  times  after  an  initial 
time  when  the  macroscopic  fields  and  polarizations  are  zero.  Nonnegative  energy  expressions  provide  a 
means  for  determining  realistic  physical  limitations,  such  as  the  upper  bounds  on  the  bandwidth  and  gain  of 
antennas  or  the  lower  bounds  on  antenna  quality  factors.  Valuable  inequalities  satisfied  by  bulk  constitutive 
parameters  as  well  as  by  the  group  and  energy  transport  velocities  in  materials  and  metamaterials  can  be 
obtained  from  nonnegative  macroscopic  energy  expressions.  Also,  knowing  the  sufficient  conditions  for  the 
validity  of  the  nonnegative  energy  expressions  lends  insight  into  the  development  and  use  of  materials  and 
metamaterials  that  may  not  be  subject  to  the  restrictions  imposed  by  these  expressions. 
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2  Energy  Supplied  to  the  Charge  Carriers 

The  energy  supplied  to  the  charge  carriers  (current  density  j)  in  a  material  volume  V  by  the  electromagnetic 
fields  (electric  field  e)  during  the  time  interval  f  -  £{j  i|  given  by 

Wje(t)-Wje(to)=  f 

j(r,  t)  ■  e(r,  t)dV dt  .  (1) 


to  v 

Assume  that  the  applied  electromagnetic  fields  are  zero  until  after  the  time  to  so  that  the  in-band  electric 
field  does  no  work  on  the  charge  carriers  until  after  to  ■  If  the  charge -current  is  generated  by  bound  charge 
carriers  whose  only  net  energy  transfer  between  them  and  their  host  material  is  in  the  form  of  heat  (out-of- 
band  energy),  then  this  work,  Wje(t )  -  Wje(to).  done  on  the  bound  charge  carriers  equals  the  change  in 
in-band  kinetic  and  potential  energy  of  the  charge  carriers  plus  any  energy  converted  into  heat.  Moreover,  if 
the  macroscopic  continuum  is  passive  in  that  there  are  no  auxiliary  sources  (such  as  capacitive  energy  stored 
in  a  compressed-spring  model  of  electric  dipoles  formed  by  equal  and  opposite  electric  charges,  or  inductive 
energy  stored  in  permanent  molecular  Amperian  magnetic  dipoles)  that  can  release  energy  upon  excitation 
by  applied  fields  to  serve  as  an  active  source  of  internal  energy  adding  to  the  in-band  energy  of  the  charge 
carriers,  then  Wje(t )  -  Wje(to)  S  0.  Choosing  the  arbitrary,  constant  initial  energy  Wje(to)  equal  to  zero, 
we  have 


Wje(t)  = 


f 


j(r,  f)  ■  e(r,  t)dV dt  >  0 


(2) 


to  v 

for  bound  charge  carriers  in  passive  material  that  can  be  lossy  as  well  as  lossless.  We  will  refer  to  a  material 
that  satisfies  the  inequality  in  (2)  for  all  V  as  being  unconditionally  passive.  The  current  j(r,  t)  in  (2)  can 
include  conduction  current  as  well  as  the  current  produced  by  the  local  motion  of  the  bound  charges  if  the 
changes  in  the  kinetic  and  potential  energy  transferred  by  the  drift  velocity  of  the  conduction  charges  are 
negligible  compared  with  the  out-of-band  heat  energy  produced  by  the  conduction  current  (as  is  normally 
the  case). 

3  Poynting’  s  Theorem  for  Dipolar  Continua 

It  can  be  proven  that  the  total  electromagnetic  power  at  time  f  entering  a  closed  surface  S  lying  in  free  space 
is  given  by  the  integral  over  S  of  the  microscopic  Poynting  vector,  namely 

-P(t)  =  -  n  ■  [e(r,  f)  *  h(r,  t)]dS  (3) 

s 

where  n  is  the  unit  normal  directed  away  from  the  volume  V  enclosed  by  S,  and  the  microscopic  magnetic 
fields  are  related  by  b  =  //oh  because  S  is  in  free-space. 

From  Maxwell’ s  equations  in  an  ideal  dipolar  continuum,  it  follows  that  the  components  of  the  electric 
and  magnetic  fields,  E  and  H,  tangential  to  an  interface  between  free  space  and  the  polarized  material  are 
continuous,  provided  there  are  no  equivalent  polarization  surface  currents  at  the  interface  of  the  polarized 
material,  that  is,  no  delta  functions  in  the  electric  and  magnetic  polarization  densities,  P  or  M,  at  the  in¬ 
terface  [5],  Such  delta  functions  can  occur  for  polarization  with  constitutive  parameters  that  are  strongly 
spatially  dispersive  [  1  ]  and  even  for  spatially  nondispersive  material  characterized  by  extreme  anisotropic 
constitutive  parameters  with  the  values  of  some  elements  approaching  zero  and  others  infinity  at  the  in- 
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terface  [5].  If  we  avoid  these  pathologically  extreme  constitutive  parameters  and  also  assume  low  enough 
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frequencies  that  spatial  dispersion  is  negligible,  then  the  tangential  components  of  the  ideal-continuum  E 
and  H  are  continuous  across  the  free-space/dipolar-material  interface.  This  implies  that  if  we  remove  an 
infinitesimally  thin  shell  of  the  ideal  continuum  containing  S  (everywhere  that  the  closed  surface  Spasses 
through  the  dipolar  continuum),  the  value  of  n  ■  (E  x  H)  will  be  continuous  across  the  resulting  free-space 
shell.  Moreover,  since  thecontinuum  fields  equal  the  microscopic  fields  in  the  free-space  shell,  thatis 

n  ■  (E  *  H)  =  n  ■  (e  x  h)  (4) 

in  the  free-space  shell,  one  finds  from  (3)  that  the  total  instantaneous  power  flowing  across  a  closed  surface 
S  in  an  ideal  dipolar  continuum  is  given  by  the  integration  of  the  continuum  Poynting  vector  P(t)  - 

—  n  [E(r,  f)  x  H(r,  t)]dS,  or  with  the  aid  of  Maxwell’s  ideal-continuum  dipolar  equations 

p(t)  =  ~  n-  [E(r,  t)  *  H(r,  t)]dS  =  g]Xr,  t)  ^  ^  +  dB(r,  t)  ^  dy 

dt  dt 

s  v 

Thus,  we  have  determined  that  in  an  ideal  continuum  the  integral  of  the  continuum  Poynting  vector 
over  a  closed  surface  S  of  a  volume  V  exactly  equals  the  instantaneous  power  flow  across  that  surface, 
provided  the  constitutive  parameters  do  not  take  on  extreme  enough  values  and  the  spatial  dispersion  is 
not  strong  enough  over  the  operational  bandwidths  to  produce  equivalent  surface  polarization  currents  at  a 
hypothetical  interface  S  between  the  polarized  material  and  free  space.  Furthermore,  the  surface  integral  of 
the  continuum  Poynting  vector  equals  the  corresponding  volume  integral  of  the  fields  on  the  right-hand  side 
of  (5). 

It  should  be  noted  that  in  some  contexts  bianisotropic  media  are  considered  to  be  spatially  dispersive 
because  the  electric  and  magnetic  fields  in  bianisotropic  media  are  proportional  to  the  curls  of  the  magnetic 
and  electric  fields,  respectively.  However,  since  this  weak  spatial  dispersion  maintains  continuous  n'(ExH) 
across  interfaces  between  the  bianisotropic  continua  and  free  space,  for  the  puiposes  of  this  work,  we  can 
include  the  weak  spatial  dispersion  of  bianisotropic  media  within  the  category  of  spatially  nondispersive 
media. 

The  foregoing  argument  that  has  led  to  Poynting’s  theorem  in  an  ideal  dipolar  continuum  and  its  inter¬ 
pretation  in  terms  of  power  flow  can  be  applied  with  a  couple  of  modifications  to  a  macroscopic  dipolar 
continuum  comprised  of  discrete  electric  and  magnetic  dipoles.  The  first  modification  is  that  the  volume  V 
with  surface  S  to  which  Poynting’s  theorem  is  applied  should  be  no  smaller  than  a  macroscopic  volume  AV 
which  is  electrically  small  but  still  contains  so  many  dipoles  that  the  averaging  of  the  fields  over  AV  gives 
to  a  good  approximation  the  continuum  fields  satisfying  Maxwell’s  equations. 

The  second  modification  involves  a  transition  layer  of  finite  thickness  at  an  interface  between  free  space 
and  the  macroscopic  dipolar  continuum  [1],  [6].  Unlike  an  ideal  dipolar  continuum,  the  tangential  E  and 
H  fields  in  a  macroscopic  dipolar  continuum  are  not  perfectly  continuous  across  an  interface  but  are  nearly 
continuous  across  a  transition  layer  of  thickness  5  containing  the  interface.  Fortunately,  under  the  condi¬ 
tions  that  the  discrete  dipolar  material  behaves  as  a  continuum,  analytical  and  numerical  results  with  discrete 
dipolar  arrays  indicate  that  the  thickness  8  of  the  transition  layer  is  on  the  order  of  the  average  separation 
distance  d  of  the  dipoles  [7],  [8].  This  means  that  the  infinitesimal  thickness  of  the  free-space  shell  about  S 
that  was  chosen  to  derive  the  above  results  for  the  Poynting  vectors  for  the  ideal  continuum  can  be  made  ap¬ 
proximately  equal  to  the  thickness  8  of  the  transition  layer  without  appreciably  changing  the  average  fields 
within  the  volume  V .  Consequently,  the  equation  (4)  holds  to  a  good  approximation  for  S  at  the  center  of 
the  6-thick  free-space  shell  that  removes  a  discrete  number  of  dipoles;  see  Fig.  1.  The  better  the  discrete 
distribution  of  dipoles  approximates  a  continuum,  the  thinner  the  transition  layer  becomes  (6  — >  0)  and 
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the  more  accurate  is  the  approximation  for  the  macroscopic  continuum.  The  only  precedent  that  we  could 
find  for  this  distinction  between  ideal  and  macroscopic  continua  was  in  Maxwell’s  original  Treatise.  This 
discovery  became  one  of  the  highlights  of  a  2014  paper  invited  to  a  special  issue  of  the  journal,  Progress  in 
Electromagnetics  Research,  commemorating  150  years  of  Maxwell’s  equations  [2],  [3]. 


Figure  1 :  Volume  V  with  its  surface  S  centered  in  a  free-space  shell  of  thickness  S  enclosing  a  large  number 
of  discrete  electric  and  magnetic  dipoles. 

The  most  important  aspect  of  (5)  is  that  both  the  surface  and  volume  integrals  in  (5)  are  equal  to  the  instan¬ 
taneous  power  flow  P(t)  through  the  closed  surface  S  into  the  volume  V,  exactly  for  the  ideal  continuum 
and  approximately  for  the  macroscopic  continuum,  respectively.  Although  free-space  shells  are  invoked  to 
prove  (5),  the  fields  in  (5)  are  the  continuum  fields  and  thus  (5)  holds  exactly  within  the  ideal  continuum,  and 
approximately  in  the  macroscopic  continuum  since  n 1  (E  x  H)  is  continuous  across  a  free-space/dipolar- 
material  interface  under  the  conditions  specified  above  and  satisfied  by  most  dipolar  metamaterial  continua 
as  well  as  natural  dipolar  material  continua. 


4  Energy  Relations  for  Dipolar  Continua 

If  the  volume  V  is  chosen  to  be  a  macroscopic  volume  AV  that  is  electrically  small  but  large  enough  to 
contain  a  great  number  of  discrete  dipoles,  then  (5)  becomes 

{  {1 

AP(f)  ~  -  n  •  [E(r,  £)  x  H(r,  t)]dS  =  ^  ■  E(r,  t ;)  +  ^  ■  H(r,  i)  dV  .  (6) 

dt  dt 

AS  AV 

We  can  also  express  the  power  AP(t)  in  terms  of  the  microscopic  fields  and  sources  within  A  V.  This  can  be 
done  by  reinstating  the  surface  AS  of  A  V  in  the  middle  of  a  free-space  shell  of  transition  layer  thickness  6 
that  removes  a  discrete  number  of  dipoles,  as  explained  in  Section  3,  where  the  volume  A  V  contains  a  large 
number  of  discrete  electric  and  magnetic  dipoles.  Since  AS  is  now  in  free  space,  the  average  macroscopic 
E  and  H  fields  there  are  equal  to  the  microscopic  e  and  h  =  b/juo  fields,  as  explained  in  Section  3. 
Thus,  combining  the  macroscopic  Poynting’s  theorem  in  (6)  with  the  corresponding  microscopic  Poynting’s 
theorem  by  means  of  the  equality  of  the  macroscopic  and  microscopic  Poynting’s  vectors  in  (4),  it  can  be 
shown  that  the  following  relationship  and  inequality  hold  for  unconditionally  passive  material 

<f1SP(r,  t)  .  SB(r,£)  '  '  1 

— —  •  E(r,  t)  -  — —  ■  M(r.  t)  AV  yedVdt 

to  01  t0AV 
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1 


(7) 


£o|eins(r,  t)  -  Eins(r,  t) |2  +  ljbins(r,  t)  -  Bins(r,  f)|2  dv  >  0 
2AV  mo 

AV 

where  the  superscripts  “ins”  refer  to  the  microscopic  and  macroscopic  fields  produced  only  by  the  micro¬ 
scopic  and  macroscopic  dipolar  sources  inside  A  V ,  respectively. 

A  dipolar  material  is  not  necessarily  unconditionally  passive  if  it  contains  permanent  dipole  moments 
whose  magnitudes  can  be  reduced  when  subject  to  applied  fields  to  create  an  effectively  active  material. 
However,  nearly  all  electric  polarization  can  be  modeled  by  either  initially  randomly  oriented  molecules 
with  permanent  electric  dipole  moments  that  stay  practically  fixed  in  magnitude  as  they  align  in  the  applied 
field,  or  by  initially  zero  electric  dipole-moment  molecules  that  distort  in  the  applied  field  to  produce  induced 
electric  dipole  moments  [9,  p.  464],  [10,  p.  162],  In  the  latter  case,  the  electric  dipoles  have  no  initial  energy 
from  which  to  draw.  In  the  former  case,  the  energy  of  formation  of  the  electric  dipoles  cannot  be  reduced 
because  the  initial  separation  distance  between  the  equal  and  opposite  charges  of  each  dipole  cannot  increase 
to  allow  their  mutual  force  to  do  negative  work  and  release  part  of  their  initial  energy  of  formation.  Thus, 
these  two  types  of  electric  dipoles  satisfy  the  criterion  of  unconditional  passivity. 

4. 1  Macroscopic  Polarization  Energy  for  Magnetic  Dipoles 

Since  magnetic  charge  does  not  exist,  magnetic  dipoles  are  created  by  molecules  with  circulating  currents 
that  can  be  modeled  by  perfectly  electrically  conducting  (PEC)  wire  loops.  If  the  wire  loops  carry  no 
permanent  current  so  that  all  the  circulating  current  and  magnetic  dipole  moments  are  induced  by  the  applied 
fields,  then  the  distribution  of  molecules  forms  a  diamagnetic  macroscopic  continuum  whose  low-frequency 
magnetic  susceptibility  is  less  than  zero. 


4.1.1  Macroscopic  Polarization  Energy  for  Diamagnetism 


In  the  case  of  diamagnetism,  the  initial  value  of  the  current  and  magnetic  dipole  moment  of  the  wire  loops 
are  zero.  Therefore,  diamagnetic -material  continua  that  can  also  contain  electric  polarization  produced 
by  electric  dipoles  with  fixed-magnitude  or  initially  zero  electric  dipole  moments  exhibit  the  property  of 
unconditional  passivity  and  thus  diamagnetic  material  satisfies  the  inequality  in  (7),  specifically 

/  1  I 

<5P(r,  f )  ,  SB(r,  t)  '  ' 


to 


di 


■  E(r,  t)  - 


dt 


■  M(r,  t ) 


dt  >  0 . 


(8) 


As  far  as  we  are  aware,  this  is  the  first  consistent,  positive  semi-definite  energy  expression  that  has  been 
derived  for  macroscopic  diamagnetic  material. 

As  a  simple  example,  consider  a  material  with  negligible  loss  and  constant  real  permittivity  E  and  perme¬ 
ability  p  over  the  operational  baseband  bandwidth  \  a\  <  ao-  Then  P  =  (E- Eo)E  and  M  =  (1/ Mo  _  1/ ,t/)B 
over  the  baseband  bandwidth,  and  (8)  reduces  to 


(E-Eo) |E|2  -  (1/mo  ~  1/^)|B|2  >  0 . 


(9) 


Since  E  and  B  can  be  chosen  equal  to  zero  independently,  this  inequality  reveals  that 


E>Eo,  0<m^Uo 

which  confirms  that  the  inequality  in  (8)  applies  to  diamagnetic  continua. 


(10) 
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5  Application  of  Energy  Relations  to  Bianisotropic  Continua 


In  this  section,  the  energy  theorems  summarized  immediately  above  are  applied  to  linear,  passive,  spatially 
nondispersive  media  to  obtain  frequency-domain  expressions  for  internal  energy  densities  in  lossless  media 
and  inequalities  that  the  linear  constitutive  relations  must  obey  in  lossless  media.  The  most  general  linear, 
spatially  nondispersive  constitutive  relations  are  those  for  bianisotropic  media  and  are  given  in  the  frequency 
domain  as 


D«(r)  =  E(r)  ■  E«(r)  +  t(i)  ■  Hw(r)  (11a) 

Boj(r)  =  m  ■  Hw(r)  +  v(r)  ■  Ew(r)  (lib) 


where  p(r),  E(r),  and  [v(r),  f(f)]  are  the  permeability  dyadic,  the  permittivity  dyadic,  and  the  magneto¬ 
electric  dyadics,  respectively.  Like  the  fields,  they  are,  in  general,  functions  of  frequency  a  and  position  r 

within  the  media. 

If  the  magnetic  polarization  of  the  bianisotropic  material  is  produced  by  diamagnetic  dipoles,  then  the 
inequality  in  (8)  applies.  It  can  be  rewritten  as 


{t1 


to 


SD(r,  f  )  .  5B(r,  t) 

■  E(r,  t )  + 


8t 


dt 


2+J0IB(r’‘)|2 


(12) 


Using  an  analysis  similar  to  the  one  in  [11]  and  [12]  but  applied  to  (12),  we  find  for  a  lossless  frequency 
window 


Re  E£,  ■  («E)'  ■  Ew  +  Hw  ■  (ojT)  '  H«  +  Ew  ■  (gj(vt  +  t  )  1  H^,  ^  £o|Ea)|2  +  ^,  ,2  (13) 


Mo 


|B«| 


which  inplies  that 


[{(oEu)  -Eo]>aEu/2>0 


(14a) 


with  a  similar  (though  not  identical)  inequality  for  the  diagonal  elements  of  the  permeability  dyadic  if  the 
orientation  of  the  xyz  coordinate  system  is  chosen  at  each  point  in  space  to  make  the  permeability  dyadic 
diagonal  (which  is  always  possible  because  the  lossless  permeability  dyadic  is  Hermitian),  namely 

'  2 

[(o/Jii)  -  jun/juo\  >  a/Ju/2  >  0  .  (14b) 


Note  that  as  a  — *  0,  (14a)  implies  that 


Eu(cd—>Q)-EQ>0 


(15) 


(assuming  lima,^o(cof?z)  =  0).  However,  (14b)  implies  that 

juu(a  —>  0)[1  -  juu(a  — ►  0)/qo]  >  0 


(16) 


or,  equivalently 


0  <  juu(od  -»  0)  <  (17) 

for  diamagnetic  material  (assuming  lirriM^oCcu^zi)  =  0).  These  inequalities  are  compatible  with  diamagnetic 
material,  yet  incompatible  with  ordinary  nondiamagnetic  (paramagnetic)  material. 


7 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Acknowledgment 

Helpful  discussions  with  Professors  M.G.  Silveirinha  and  A.  Alu  are  gratefully  acknowledged. 

References 

[1]  A.D.  Yaghjian,  A.  Alu,  and  M.G.  Silveirinha.  Anisotropic  representation  for  spatially  dispersive  peri¬ 

odic  metamaterial  arrays,  ch.  13  in  Transformation  Electromagnetics  and  Metamaterials,  D.H.  Werner, 
and  D.-H.  Kwon,  Eds.,  Springer,  New  York,  2014;  see  also _ .  Homogenization  of  spatially  dis¬ 

persive  metamaterial  arrays  in  terms  of  generalized  electric  and  magnetic  polarizations.  Photonics  and 
Nanostructures,  1 1, 374-396, 2013. 

[2]  A.D.  Yaghjian,  “Reflections  on  Maxwell’s  Treatise,”  PIER ,  vol.  149,  pp.  217-249,  November  2014. 

[3]  A.D.  Yaghjian,  “An  overview  of  Maxwell’s  Treatise,”  FERMAT,  vol.  11,  September-October,  2015. 

[4]  A.D.  Yaghjian  “Power  &  energy  relations  for  macroscopic  dipolar  continua,”  IEEE-APS/URSI  Sym¬ 
posium  Digest ,  p.  367,  July  2015. 

[5]  A.D.  Yaghjian,  “Extreme  electromagnetic  boundary  conditions  and  their  manifestation  at  the  inner  sur¬ 
faces  of  spherical  and  cylindrical  cloaks,”  Metamaterials  Journal,  vol.  4,  pp.  70-76,  August-September 
2010. 

[6]  A.D.  Yaghjian,  “Boundary  conditions  for  electric  quadrupolar  continua,”  Radio  Science,  vol.  49,  pp. 
1289-1299,  December  2014. 

[7]  C.R.  Simovski  and  S.A.  Tretyakov,  “Local  constitutive  parameters  of  metamaterials  from  an  effective- 
medium  perspective,”  Phys.  Rev.  B,  vol.  75,  pp.  19511 1(1-10),  2007. 

[8]  A.D.  Scher  and  E.F.  Kuester,  “Boundary  effects  in  the  electromagnetic  response  of  a  metamaterial  in 
the  case  of  normal  incidence,”  PIER  B,  vol.  14,  pp.  341-381, 2009. 

[9]  C.  Kittel,  Introduction  to  Solid  State  Physics,  8th  Ed. ,  Wiley,  2005. 

[10]  J.D.  Jackson,  Classical  Electrodynamics,  3rd  Edition,  New  York:  Wiley,  1999. 

[11]  A.D.  Yaghjian  and  S.R.  Best,  “Impedance,  bandwidth,  and  Q  of  antennas,”  IEEE  Trcms.  Antennas 
Propagat.,  vol.  53,  pp.  1298-1324,  April  2005;  Correction,  vol.  55,  p.  3748,  December  2007. 

[12]  A.D.  Yaghjian,  “Internal  energy,  Q-energy,  Poynting’s  theorem,  and  the  stress  dyadic  in  dispersive 
material,”  IEEE  Trans.  Antennas  Propagat.,  vol.  55,  pp.  1495-1505,  June  2007. 


8 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other 
aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other 
provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

01/06/2016 


2.  REPORT  TYPE 

Final 


3.  DATES  COVERED  (From  -  To) 

15  March  2013-14  March  2016 


4.  TITLE  AND  SUBTITLE 

Energy  Relations  in  Natural  and  Artificial  Diamagnetic  Materials 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

FA9550-1 3-1 -0033 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Arthur  D.  Yaghjian 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

S4  Inc. 

209  Burlington  Road,  Suite  105 
Bedford,  MA01730 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Office  of  Scientific  Research 
875  North  Randolph  Street,  RM  3112 
Arlington,  VA  22203 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL/AFOSR/RSE 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Distribution  Unlimited 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

Positive  semi-definite  (nonnegative)  expressions  for  the  time-domain  macroscopic  energy  density  in  passive,  spatially 
nondispersive,  diamagnetic,  dipolar  continua  are  derived  from  the  underlying  microscopic  Maxwellian  equations  satisfied  by 
classical  models  of  discrete  bound  diamagnetic  molecules  or  inclusions  of  the  material  or  metamaterial  continua.  The 
macroscopic  diamagnetic  energy  expressions  are  applied  to  linear,  passive,  spatially  nondispersive  diamagnetic  media  to 
obtain  frequency-domain  expressions  for  internal  energy  densities  in  lossless  diamagnetic  media  and  inequalities  that  the 
linear  constitutive  relations  must  obey  in  lossless  diamagnetic  media. 


15.  SUBJECT  TERMS 

Energy  Relations,  Diamagnetic  Materials  and  Metamaterials,  Macroscopic  and  Microscopic  Maxwellian  Equations 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

Arje  Nachman 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

u 

u 

u 

SAR 

7 

703-696-8427 

DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


Final  Performance  Report  for  Grant  FA9550-13-1-0033: 
“Energy  Relations  in  Natural  and  Artificial  Diamagnetic  Materials” 

Program  Manager:  Dr.  Arje  Nachman,  Principal  Investigator:  Dr.  Arthur  D.  Yaghjian 


1  Introduction 

As  a  consequence  of  the  insight  gained  from  the  rigorous  homogenization  theory  developed  by  the  principal 
investigator  in  collaboration  with  Andrea  Alu  and  Mario  Sil  veirinha  [  1  ]  for  the  magnetic  permeability  as  well 
as  electric  permittivity  of  3D  spatially  dispersive  periodic  arrays,  the  AFOSR  three-year  Grant  FA9550-13- 
1-0033  for  $1 10,889  and  its  Additional  Work  Effort  for  $72,600  was  awarded  to  S4-Inc  to  have  Dr.  Yaghjian 
apply  this  theory  to  determine  consistent  expressions  for  the  power  and  energy  densities  in  a  diamagnetic 
continuum.  The  tasks  proposed  and  successfully  completed  for  the  three  years  were  in  brief: 

1 .  To  evaluate  the  external  and  internal  power  and  energy  relations  for  a  realistic  model  of  a  diamagnetic 
molecule  or  inclusion,  such  as  a  perfectly  electrically  conducting  (PEC)  wire  loop; 

2.  To  determine  a  reliable  method  for  obtaining  macroscopic  continuum  relations  of  power  and  en¬ 
ergy  from  the  microscopic  Maxwellian  equations  describing  the  power  and  energy  of  diamagnetic 
molecules  or  inclusions. 

3.  To  incorporate  the  results  for  the  macroscopic  power  and  energy  relations  into  the  continuum  formu¬ 
lation  for  diamagnetic  materials  or  metamaterials; 

4.  To  combine  the  diamagnetic  relations  for  power  and  energy  with  those  of  electric  polarization  to 
obtain  positive  semi-definite  (that  is,  nonnegative)  energy  expressions  for  diamagnetic  continua,  then 
document  the  results  of  the  research  in  archival  publications. 

These  accomplishments,  which  were  recently  submitted  to  the  journal  Physical  Review  B,  will  be  summa¬ 
rized  in  the  remainder  of  this  final  report.  Some  of  the  results  can  also  be  found  in  the  references  [2] — [4] . 

There  are  several  reasons  for  determining  nonnegative  macroscopic  energies  from  the  microscopic 
Maxwellian  equations  for  classical  diamagnetic  dipolar  continua.  First,  it  is  one  of  the  remaining  unre¬ 
solved  problems  in  classical  electromagnetic  theory  even  though  the  problem  is  relatively  easy  to  state: 
given  a  volume  of  a  macroscopic  continuum  satisfying  Maxwell’s  dipolar  equations  and  illuminated  by  ex¬ 
ternal  fields,  are  there  macroscopic  energy  densities  that  never  become  negative  for  all  times  after  an  initial 
time  when  the  macroscopic  fields  and  polarizations  are  zero.  Nonnegative  energy  expressions  provide  a 
means  for  determining  realistic  physical  limitations,  such  as  the  upper  bounds  on  the  bandwidth  and  gain  of 
antennas  or  the  lower  bounds  on  antenna  quality  factors.  Valuable  inequalities  satisfied  by  bulk  constitutive 
parameters  as  well  as  by  the  group  and  energy  transport  velocities  in  materials  and  metamaterials  can  be 
obtained  from  nonnegative  macroscopic  energy  expressions.  Also,  knowing  the  sufficient  conditions  for  the 
validity  of  the  nonnegative  energy  expressions  lends  insight  into  the  development  and  use  of  materials  and 
metamaterials  that  may  not  be  subject  to  the  restrictions  imposed  by  these  expressions. 
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2  Energy  Supplied  to  the  Charge  Carriers 


The  energy  supplied  to  the  charge  carriers  (current  density  j)  in  a  material  volume  V  by  the  electromagnetic 
fields  (electric  held  e)  during  the  time  interval  t  —  t o  is  given  by 

t 

Wje(t)  -  Wje(to)  =  J  J  j(r,  t')  •  e(r,  t')dVdt' .  (1) 

t0  v 

Assume  that  the  applied  electromagnetic  fields  are  zero  until  after  the  time  to  so  that  the  in-band  electric 
held  does  no  work  on  the  charge  carriers  until  after  to ■  If  the  charge-current  is  generated  by  bound  charge 
carriers  whose  only  net  energy  transfer  between  them  and  their  host  material  is  in  the  form  of  heat  (out-of- 
band  energy),  then  this  work,  Wje{t)  —  Wje(to),  done  on  the  bound  charge  carriers  equals  the  change  in 
in-band  kinetic  and  potential  energy  of  the  charge  carriers  plus  any  energy  converted  into  heat.  Moreover,  if 
the  macroscopic  continuum  is  passive  in  that  there  are  no  auxiliary  sources  (such  as  capacitive  energy  stored 
in  a  compressed-spring  model  of  electric  dipoles  formed  by  equal  and  opposite  electric  charges,  or  inductive 
energy  stored  in  permanent  molecular-  Amperian  magnetic  dipoles)  that  can  release  energy  upon  excitation 
by  applied  helds  to  serve  as  an  active  source  of  internal  energy  adding  to  the  in-band  energy  of  the  charge 
carriers,  then  Wje(t)  —  Wje(to)  >  0.  Choosing  the  arbitrary,  constant  initial  energy  Wje(to)  equal  to  zero, 
we  have 

t 

Wje(t )  =  j  J  j  (r ,  t')  ■  e(r,  t')dVdt'  >  0  (2) 

t0  V 

for  bound  charge  carriers  in  passive  material  that  can  be  lossy  as  well  as  lossless.  We  will  refer  to  a  material 
that  satisfies  the  inequality  in  (2)  for  all  V  as  being  unconditionally  passive.  The  current  j(r,  t)  in  (2)  can 
include  conduction  current  as  well  as  the  current  produced  by  the  local  motion  of  the  bound  charges  if  the 
changes  in  the  kinetic  and  potential  energy  transferred  by  the  drift  velocity  of  the  conduction  charges  are 
negligible  compared  with  the  out-of-band  heat  energy  produced  by  the  conduction  current  (as  is  normally 
the  case). 

3  Poynting’s  Theorem  for  Dipolar  Continua 

It  can  be  proven  that  the  total  electromagnetic  power  at  time  t  entering  a  closed  surface  S  lying  in  free  space 
is  given  by  the  integral  over  S  of  the  microscopic  Poynting  vector,  namely 

P(t)  =  —  J  n  •  [e(r,  t )  x  h(r,  t)]dS  (3) 

s 

where  n  is  the  unit  normal  directed  away  from  the  volume  V  enclosed  by  S,  and  the  microscopic  magnetic 
fields  are  related  by  b  =  //qIi  because  S  is  in  free-space. 

From  Maxwell’s  equations  in  an  ideal  dipolar  continuum,  it  follows  that  the  components  of  the  electric 
and  magnetic  fields,  E  and  H,  tangential  to  an  interface  between  free  space  and  the  polarized  material  are 
continuous,  provided  there  are  no  equivalent  polarization  surface  currents  at  the  interface  of  the  polarized 
material,  that  is,  no  delta  functions  in  the  electric  and  magnetic  polarization  densities,  P  or  M,  at  the  in¬ 
terface  [5].  Such  delta  functions  can  occur  for  polarization  with  constitutive  parameters  that  are  strongly 
spatially  dispersive  [1]  and  even  for  spatially  nondispersive  material  characterized  by  extreme  anisotropic 
constitutive  parameters  with  the  values  of  some  elements  approaching  zero  and  others  infinity  at  the  in¬ 
terface  [5].  If  we  avoid  these  pathologically  extreme  constitutive  parameters  and  also  assume  low  enough 
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frequencies  that  spatial  dispersion  is  negligible,  then  the  tangential  components  of  the  ideal-continuum  E 
and  H  arc  continuous  across  the  free-space/dipolar-material  interface.  This  implies  that  if  we  remove  an 
infinitesimally  thin  shell  of  the  ideal  continuum  containing  S  (everywhere  that  the  closed  surface  S  passes 
through  the  dipolar  continuum),  the  value  of  n  •  (E  x  H)  will  be  continuous  across  the  resulting  free-space 
shell.  Moreover,  since  the  continuum  fields  equal  the  microscopic  fields  in  the  free-space  shell,  that  is 

n  •  (E  x  H)  =  n  •  (e  x  h)  (4) 


in  the  free-space  shell,  one  finds  from  (3)  that  the  total  instantaneous  power  flowing  across  a  closed  surface 
S  in  an  ideal  dipolar  continuum  is  given  by  the  integration  of  the  continuum  Poynting  vector  P(t)  = 
~  fsA  '  [E(r> l)  x  H(r,  t)]dS,  or  with  the  aid  of  Maxwell's  ideal-continuum  dipolar  equations 


n  •  [E(r,  t)  x  H(r,  t)]dS 


v 


<9D(r,  t) 

Ft 


d,V. 


(5) 


Thus,  we  have  determined  that  in  an  ideal  continuum  the  integral  of  the  continuum  Poynting  vector 
over  a  closed  surface  S  of  a  volume  V  exactly  equals  the  instantaneous  power  flow  across  that  surface, 
provided  the  constitutive  parameters  do  not  take  on  extreme  enough  values  and  the  spatial  dispersion  is 
not  strong  enough  over  the  operational  bandwidths  to  produce  equivalent  surface  polarization  currents  at  a 
hypothetical  interface  S  between  the  polarized  material  and  free  space.  Furthermore,  the  surface  integral  of 
the  continuum  Poynting  vector  equals  the  corresponding  volume  integral  of  the  fields  on  the  right-hand  side 
of  (5). 

It  should  be  noted  that  in  some  contexts  bianisotropic  media  are  considered  to  be  spatially  dispersive 
because  the  electric  and  magnetic  fields  in  bianisotropic  media  are  proportional  to  the  curls  of  the  magnetic 
and  electric  fields,  respectively.  However,  since  this  weak  spatial  dispersion  maintains  continuous  n-(ExH) 
across  interfaces  between  the  bianisotropic  continua  and  free  space,  for  the  purposes  of  this  work,  we  can 
include  the  weak  spatial  dispersion  of  bianisotropic  media  within  the  category  of  spatially  nondispersive 
media. 

The  foregoing  argument  that  has  led  to  Poynting’s  theorem  in  an  ideal  dipolar  continuum  and  its  inter¬ 
pretation  in  terms  of  power  flow  can  be  applied  with  a  couple  of  modifications  to  a  macroscopic  dipolar 
continuum  comprised  of  discrete  electric  and  magnetic  dipoles.  The  first  modification  is  that  the  volume  V 
with  surface  S  to  which  Poynting’s  theorem  is  applied  should  be  no  smaller  than  a  macroscopic  volume  AV 
which  is  electrically  small  but  still  contains  so  many  dipoles  that  the  averaging  of  the  fields  over  AV  gives 
to  a  good  approximation  the  continuum  fields  satisfying  Maxwell’s  equations. 

The  second  modification  involves  a  transition  layer  of  finite  thickness  at  an  interface  between  free  space 
and  the  macroscopic  dipolar  continuum  [1],  [6].  Unlike  an  ideal  dipolar  continuum,  the  tangential  E  and 
H  fields  in  a  macroscopic  dipolar  continuum  are  not  perfectly  continuous  across  an  interface  but  are  nearly 
continuous  across  a  transition  layer  of  thickness  S  containing  the  interface.  Fortunately,  under  the  condi¬ 
tions  that  the  discrete  dipolar  material  behaves  as  a  continuum,  analytical  and  numerical  results  with  discrete 
dipolar  arrays  indicate  that  the  thickness  5  of  the  transition  layer  is  on  the  order  of  the  average  separation 
distance  d  of  the  dipoles  [7],  [8].  This  means  that  the  infinitesimal  thickness  of  the  free-space  shell  about  S 
that  was  chosen  to  derive  the  above  results  for  the  Poynting  vectors  for  the  ideal  continuum  can  be  made  ap¬ 
proximately  equal  to  the  thickness  5  of  the  transition  layer  without  appreciably  changing  the  average  fields 
within  the  volume  V.  Consequently,  the  equation  (4)  holds  to  a  good  approximation  for  S  at  the  center  of 
the  (5-thick  free-space  shell  that  removes  a  discrete  number  of  dipoles;  see  Fig.  1.  The  better  the  discrete 
distribution  of  dipoles  approximates  a  continuum,  the  thinner  the  transition  layer  becomes  (5  — >  0)  and 
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the  more  accurate  is  the  approximation  for  the  macroscopic  continuum.  The  only  precedent  that  we  could 
find  for  this  distinction  between  ideal  and  macroscopic  continua  was  in  Maxwell’s  original  Treatise.  This 
discovery  became  one  of  the  highlights  of  a  2014  paper  invited  to  a  special  issue  of  the  journal.  Progress  in 
Electromagnetics  Research,  commemorating  150  years  of  Maxwell’s  equations  [2],  [3]. 


Figure  1:  Volume  V  with  its  surface  S  centered  in  a  free-space  shell  of  thickness  5  enclosing  a  large  number 
of  discrete  electric  and  magnetic  dipoles. 

The  most  important  aspect  of  (5)  is  that  both  the  surface  and  volume  integrals  in  (5)  are  equal  to  the  instan¬ 
taneous  power  flow  P(t)  through  the  closed  surface  S  into  the  volume  V,  exactly  for  the  ideal  continuum 
and  approximately  for  the  macroscopic  continuum,  respectively.  Although  free-space  shells  are  invoked  to 
prove  (5),  the  fields  in  (5)  are  the  continuum  fields  and  thus  (5)  holds  exactly  within  the  ideal  continuum,  and 
approximately  in  the  macroscopic  continuum  since  n  •  (E  x  H)  is  continuous  across  a  free-space/dipolar- 
material  interface  under  the  conditions  specified  above  and  satisfied  by  most  dipolar  metamaterial  continua 
as  well  as  natural  dipolar  material  continua. 


4  Energy  Relations  for  Dipolar  Continua 


If  the  volume  V  is  chosen  to  be  a  macroscopic  volume  AV  that  is  electrically  small  but  large  enough  to 
contain  a  great  number  of  discrete  dipoles,  then  (5)  becomes 


A P(t)  «  —  J  n  •  [E(r,  t)  x  H(r,  t)]dS  =  j 

AS  AV 


<9D(r,  t) 
dt 


E(M)  +  ^A.H(M) 


dV .  (6) 


We  can  also  express  the  power  AP(t)  in  terms  of  the  microscopic  fields  and  sources  within  AV.  This  can  be 
done  by  reinstating  the  surface  AS  of  A V  in  the  middle  of  a  free-space  shell  of  transition  layer  thickness  5 
that  removes  a  discrete  number  of  dipoles,  as  explained  in  Section  3,  where  the  volume  AV  contains  a  large 
number  of  discrete  electric  and  magnetic  dipoles.  Since  A  .S'  is  now  in  free  space,  the  average  macroscopic 
E  and  H  fields  there  are  equal  to  the  microscopic  e  and  h  =  h/po  fields,  as  explained  in  Section  3. 
Thus,  combining  the  macroscopic  Poynting’s  theorem  in  (6)  with  the  corresponding  microscopic  Poynting’s 
theorem  by  means  of  the  equality  of  the  macroscopic  and  microscopic  Poynting’s  vectors  in  (4),  it  can  be 
shown  that  the  following  relationship  and  inequality  hold  for  unconditionally  passive  material 


/ 


<9P(r,  t') 

dt’ 


E(r,  t') 


<9B(r,  t') 

dt’ 


M(r,  t') 


dt' 


t. 


t0AV 


•  edV  dt' 
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/U|eins(r,t)  -  Eins(r,  t)\2  +  -|bins(r,  t)  -  Bins(r,f)|2l  dV  >  0  (7) 

AIM/  J  l  Ho 

av 

where  the  superscripts  “ins”  refer  to  the  microscopic  and  macroscopic  fields  produced  only  by  the  micro¬ 
scopic  and  macroscopic  dipolar  sources  inside  AV,  respectively. 

A  dipolar  material  is  not  necessarily  unconditionally  passive  if  it  contains  permanent  dipole  moments 
whose  magnitudes  can  be  reduced  when  subject  to  applied  fields  to  create  an  effectively  active  material. 
However,  nearly  all  electric  polarization  can  be  modeled  by  either  initially  randomly  oriented  molecules 
with  permanent  electric  dipole  moments  that  stay  practically  fixed  in  magnitude  as  they  align  in  the  applied 
field,  or  by  initially  zero  electric  dipole-moment  molecules  that  distort  in  the  applied  field  to  produce  induced 
electric  dipole  moments  [9,  p.  464],  [10,  p.  162].  In  the  latter  case,  the  electric  dipoles  have  no  initial  energy 
from  which  to  draw.  In  the  former  case,  the  energy  of  formation  of  the  electric  dipoles  cannot  be  reduced 
because  the  initial  separation  distance  between  the  equal  and  opposite  charges  of  each  dipole  cannot  increase 
to  allow  their  mutual  force  to  do  negative  work  and  release  part  of  their  initial  energy  of  formation.  Thus, 
these  two  types  of  electric  dipoles  satisfy  the  criterion  of  unconditional  passivity. 

4.1  Macroscopic  Polarization  Energy  for  Magnetic  Dipoles 

Since  magnetic  charge  does  not  exist,  magnetic  dipoles  are  created  by  molecules  with  circulating  currents 
that  can  be  modeled  by  perfectly  electrically  conducting  (PEC)  wire  loops.  If  the  wire  loops  carry  no 
permanent  current  so  that  all  the  circulating  current  and  magnetic  dipole  moments  are  induced  by  the  applied 
fields,  then  the  distribution  of  molecules  forms  a  diamagnetic  macroscopic  continuum  whose  low-frequency 
magnetic  susceptibility  is  less  than  zero. 


4.1.1  Macroscopic  Polarization  Energy  for  Diamagnetism 

In  the  case  of  diamagnetism,  the  initial  value  of  the  current  and  magnetic  dipole  moment  of  the  wire  loops 
are  zero.  Therefore,  diamagnetic-material  continua  that  can  also  contain  electric  polarization  produced 
by  electric  dipoles  with  fixed-magnitude  or  initially  zero  electric  dipole  moments  exhibit  the  property  of 
unconditional  passivity  and  thus  diamagnetic  material  satisfies  the  inequality  in  (7),  specifically 


FlW,,  9B(r,t') 
qP  ■  E(r,  t )  — 


M(r,  t')  dt'  >  0  . 


As  far  as  we  are  aware,  this  is  the  first  consistent,  positive  semi-definite  energy  expression  that  has  been 
derived  for  macroscopic  diamagnetic  material. 

As  a  simple  example,  consider  a  material  with  negligible  loss  and  constant  real  permittivity  e  and  perme¬ 
ability  p  over  the  operational  baseband  bandwidth  |oz|  <  cuo-  Then  P  =  (e  —  eo)E  and  M  =  ( 1  / /zo  —  l//x)B 
over  the  baseband  bandwidth,  and  (8)  reduces  to 

(e-e0)|E|2-(l//i0-l//r)|B|2>0.  (9) 

Since  E  and  B  can  be  chosen  equal  to  zero  independently,  this  inequality  reveals  that 

e  >  eo  ,  0  <  p  <  po  (10) 

which  confirms  that  the  inequality  in  (8)  applies  to  diamagnetic  continua. 
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5  Application  of  Energy  Relations  to  Bianisotropic  Continua 


In  this  section,  the  energy  theorems  summarized  immediately  above  are  applied  to  linear,  passive,  spatially 
nondispersive  media  to  obtain  frequency-domain  expressions  for  internal  energy  densities  in  lossless  media 
and  inequalities  that  the  linear  constitutive  relations  must  obey  in  lossless  media.  The  most  general  linear, 
spatially  nondispersive  constitutive  relations  are  those  for  bianisotropic  media  and  are  given  in  the  frequency 
domain  as 

Du,(r)  =  e(r)  •  Ew(r)  +  r(r)  •  H^(r)  (11a) 

B^(r)  =  7Z(r)  •  H^r)  +F(r)  •  E^(r)  (lib) 

where  /l(r),  e(r),  and  [I7(r),  r(r)]  are  the  permeability  dyadic,  the  permittivity  dyadic,  and  the  magneto¬ 
electric  dyadics,  respectively.  Like  the  fields,  they  are,  in  general,  functions  of  frequency  u  and  position  r 
within  the  media. 

If  the  magnetic  polarization  of  the  bianisotropic  material  is  produced  by  diamagnetic  dipoles,  then  the 
inequality  in  (8)  applies.  It  can  be  rewritten  as 

t 


to 

Using  an  analysis  similar  to  the  one  in  [11]  and  [12]  but  applied  to  (12),  we  find  for  a  lossless  frequency 
window 

Re  |e*  •  (cue)'  ■  Ew  +  H*  •  (up)'  •  +  Ew  •  [(lu(vt  +  r*)]'  •  H^}  > 

which  inplies  that 

[(coeu)'  -  e0]  >  coe'u/2  >  0  (14a) 

with  a  similar  (though  not  identical)  inequality  for  the  diagonal  elements  of  the  permeability  dyadic  if  the 
orientation  of  the  xyz  coordinate  system  is  chosen  at  each  point  in  space  to  make  the  permeability  dyadic 
diagonal  (which  is  always  possible  because  the  lossless  permeability  dyadic  is  Hermitian),  namely 

[(uj/MiY  -  Vu/v 0]  >  un'u/ 2  >  0  .  (14b) 

Note  that  as  u  — >  0,  (14a)  implies  that 


eo  I  |  H - 1  Bu 

Ato 


(13) 


8D(M)-E(M')+8B<r'< 


dtJ 


dt' 


H(r,  t') 


dt’>- 


e0|E(r,f)|2  +  — jB(r,f)|2 
Vo 


(12) 


ea((v  ->  0)  -  e0  >  0 

(assuming  lim^^oC^e^)  =  0).  However,  (14b)  implies  that 

yu(LO  ->  0) [1  -  mi(uj  -»•  0 )/n0\  >  0 


(15) 


(16) 


or,  equivalently 


0  <  Vlliu  ->•  0)  <  Vo 


(17) 


for  diamagnetic  material  (assuming  lirn^^o  (ljJVu)  =  0)-  These  inequalities  are  compatible  with  diamagnetic 
material,  yet  incompatible  with  ordinary  nondiamagnetic  (paramagnetic)  material. 
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